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O ■ Abstract 

o ■ 

' Oscillations of active to sterile neutrinos with a small mixing angle sin 26 < 

r"| ' 10"'^ could generate a large lepton asymmetry in the Early Universe. The final 

Q-i! order of magnitude of the lepton asymmetry rj is mainly determined by its 

Q^' growth in the last stage of evolution, the so called power-law regime. There 

^ ■ exist two contradictory results in the literature, rj oc and t] oc T~^, where T 

• • , is the background medium temperature. In the first case, the lepton asymmetry 

I does not exceed values of 10~^ for \6m'^\ < leV^, while in the second case it 

^ ■ can become larger than 10~^. In this work we analytically investigate the case 

rj oc T^^, using a new approach to solve the kinetic equations. We find that 
the power-law solution rj oc is not self-consistent. Instead, we find the 
power law r] oc T~^^^^ to be a good approximation, which leads to a large final 
asymmetry. 
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1 Introduction 



Neutrinos play an important role during Big Bang Nucleosynthesis (BBN) due to 
their contribution to the total energy density. Furthermore, the electron neutrinos 
directly influence the neutron-proton reactions. In this context, the Standard Model 
with its three active massless neutrinos agrees very well with the observations of the 
primordial light element abundances: a recent analysis ^ of observational data claims 
that any particles beyond the Standard Model must not contribute more than 20% 
to the total energy density of a single neutrino flavor. 

The recent discovery of u^-Ur oscillations at the Super-Kamiokande detector 
has proven the existence of neutrino masses. Reasonable neutrino masses of order 
nil, < 1 eV will not spoil the standard BBN picture, because BBN takes place at 
temperatures T ~ 0.1 MeV ^ m,y. Also, mixing between the active neutrino flavors 
plays no role since they are equally populated during the relevant epoch. 

But the situation dramatically changes if one supposes the existence of one or more 
sterile neutrinos: even a partial excitation of more than 20% of a sterile neutrino state, 
e.g. via oscillations, will contradict the observational data of light element abundances. 
Actually, introducing sterile neutrinos is quite natural, since right-handed neutrinos 
would necessarily be sterile in the Standard Model sense! 

Another important place in which sterile neutrinos can play a role is the creation 
of lepton asymmetry in the active neutrino sector. For example, if a significant 
asymmetry between z/e and is established, this could change the neutron-proton 
ratio and lead to a contradiction with the observational data. 

Early considerations of the simplified kinetic equations found no significant in- 
crease of the asymmetry in the case where the sterile neutrino has a mass m,j^ <^ 
1 MeV and any vacuum mixing angle and concluded that this asymmetry was always 
small, ?7 ^ 1 0, However, later it was found that in the region of small mixing 
angles the asymmetry can increase significantly and reach values (9(0.1) 0. This 
statement was confirmed by numerical calculations [^] where a large final value of the 
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asymmetry was found for the first time. Recently, more accurate numerical calcu- 
lations also confirmed the increase of asymmetry to large values (9(0.1) [0, §]. The 
reason for such a big increase of the asymmetry is that it grows according to a power 
law close to T~^. 

Contrary to this statement is the result of an analytical approach by Dolgov et al. 
IP], who found an additional counterterm in the effective equation for the asymmetry, 
such that the power law would change from to T~^. Thus the final asymmetry 
would be 0(10"^). 

In this paper, we analyze the method of and show that it is valid only when 
the evolution of the kinetic equations dictated by the collision terms, i.e. when the 
temperature of the primordial plasma is large enough, T ^ 1 MeV. For smaller 
temperatures, T < 2-3 MeV, collisions become less important, and the evolution 
of the neutrinos can be described by the well-known Mikheyev-Smirnov-Wolfenstein 
(MSW) effect [|l^. In this small temperature region the approach of is not valid. 
This effect has been studied in 0. 

The counterterm found in which leads to the power law T~^, only becomes 
compatible to the main term when the applied perturbative expansion breaks down. 
This fact suggests that the power law T"^ could be wrong, but strictly speaking it does 
not prove this. Therefore, we develop a new independent analytical method which can 
be applied in the regions where the method of does not work. Instead of expanding 
in the coefficients of the kinetic equations, we solve these equations pertubatively for 
the case when the creation rate of sterile neutrinos is not significantly large. Using 
this new method we show that the power law is not a self-consistent solution 
for the lepton asymmetry. Instead, the power law T^^^/'^ is approximately valid, and 
therefore the asymmetry can reach large values of order 0.1, although our new method 
in its present form does not allow us to calculate the final value of the asymmetry. 

In section 2 we define the kinetic equations. Then in section 3 we shortly discuss 
the evolution of the lepton asymmetry before it reaches the power-law regime and 
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explain the power-law solution known in literature. In section 4 we discuss some 
recent papers on this topic, especially the one claiming the power law 0. In 
section 5 we present an improved approach similar to and discuss its validity. In 
section 6 we present a new analytic solution of the kinetic equations which we use to 
check the validity of the T^^ power law. In section 7 we conclude and summarize our 
results. In the Appendix we comment on some numerical fits. 

2 Kinetic equations 

In order to describe the evolution of the system of oscillating active u^, a = e, /i, r, 
and sterile neutrinos z/g with a mass difference and a vacuum mixing angle 6 in 



the Early Universe we need to consider the density matrix formalism |11, 12] and take 



into account second order effects in the Fermi coupling constant |TT], O, KM. Then 



the following set of equations for the components of the density matrix can be found: 

i{dt- Hpdp)paa = Foipsa- Pas)/'^-i^oiPaa- feq) , 

i{dt - Hpdp)pss = -Fo{psa - Pas)/2 , 

i{dt - Hpdp)pas = WoPas + Fo{pss - paa)/'^ " i^lPas , 

i{dt - Hpdp)psa = -WoPsa - Fo{pss - Paa)/'2 - iTiPsa , (1) 

where a and s mean "active" and "sterile" respectively, Fq = dm"^ sm26/2E, Wq = 
5m^ cos26'/2E + V^^, H = .JSTiptot/Sml^ is the Hubble parameter, mpi = 1.22 x 10^^ 
MeV is the Planck mass. Furthermore, since the mass of the neutrinos is negligible in 
comparison with their momentum, the neutrino momentum and energy equal, p = E. 
Finally, feq is the equilibrium Fermi distribution function, given by 

where T is the photon temperature and p is the chemical potential for the active 
neutrinos. 
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The effective potential for (anti)-neutrinos is [15 



= TC.vGpT' + , (3) 

where E is the neutrino energy, T is the temperature of the plasma, Gp = 1.166 x 10~^ 
GeV~^ is the Fermi coupling constant, a = 1/137 is the fine structure constant, 
Ci = 0.345, and C2 = 0.61 for Ue-i^s mixing and O^"'^ = 0.17 for v^^t-i-'s mixing. 
The signs "=f" refer to neutrinos and anti-neutrinos, respectively. The individual 
contributions to the effective asymmetry rj from different particle species are 

T] = 27]^^ + 7]^^ + 7]^^ +r]e- ^n/2 for U^, 

V = 2r7^,, + r]^^ + T]y^ - r/n/2 for z/^, (4) 

and 7] for Vr is obtained from Eq. by the interchange fi ^ t. 

Equations (|l]) account exactly for the first-order terms described by the refractive 
index, while the second order terms describing the coherence breaking are approxi- 



mately modelled by the damping coefficients Tj. The latter are equal to |T6 



To = 2ri = g^(p)^^^GlT'p [1 T ZaVu. + Oivlj] , (5) 

where the expression in brackets stems from the chemical potential, and again "=p" 
denotes u and u, respectively. The coefficients of the expansion in powers of the 



asymmetry are small, = 0.1 and z^^^ = 0.04 [0, and in the following we will 
neglect these terms. The coefficients ga{p) are in general momentum-dependent, but 
in order to simplify Eqs. (||) we can take their momentum averaged values, which are 
g^^ = 3.56 and g^^ ^ = 2.5 [|. 

The anti- neutrino density matrix satisfies a set of equations similar to Eqs. (|l|), 
but with opposite sign in the antisymmetric term of and with a slight difference 
in the damping factors that is proportional to the lepton asymmetry. 

For convenience, we will introduce new variables, following the notation of Dolgov 
et al. to a certain degree, and we will restrict ourselves to treating only. First 
we note that for 1 MeV < T < 100 MeV the universe is radiation dominated so that 
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the scale factor evolves as a{t) cx T~^. This implies H = ~T/T. Therefore, if we use 
the comoving frame given by the variables x = IMeV/T and y = p/T, the LHS of 
Eqs. (0) simplify to ixHdxPij- Using ptot = ^9*T^ with the constant = 10.75, we 
get H ~ 5.5r^/mpi. 

Furthermore, we substitute 

Paa = fcq{y)[l + a{x,y)] , 
Pss = /eq(l/)[l + S{x,y)] , 

Pas=P*sa = fcq{y)[h{x,y) + il{x,y)] . (6) 

Next, we introduce a new time variable q which is defined so that for rj = the 
resonance condition, given by Wq = 0, is fulfilled when q = y. We find 



\Sm'^ \ cos 26 



x^ ^ 6.6 X lOy \6mly\ cos 20 x^ , (7) 



\\ 2C|G^a-i MeV^ 

where Smly = 5m? / eV^. For the evolution of a single momentum mode, we will use 
the variable t = q/y. 

With all these modifications, the Eqs. (|I]) have the form 

s' = Fl , 

a' = -Fl - 27a , 

h' = Wl--fh , 

I' = ^(a-s)-Wh-^l , (8) 

where the prime means differentiation with respect to r. The coefficients are derived 
using W = Wo \ 7 = Fi (^Hx^^ \ and so forth. We also split W into its 

symmetric and anti-symmetric parts, so W = U ^ Vrj. Then 

F = -Q tan 26 ^ -Q sin 26, 7 = 5/r^ U = Q{1/t^ - 1), V = Voy'^^^r-^^^ . 

(9) 
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The dimensionless constants are given by 




5.6 X 10^i/|5m2|cos2^ , 



(10) 



(11) 



mpiCiGf MeV f \6m^\ cos 26 



) 



1/6 



5.5x3 V2C|G|a-iMeV' 



5.6 X 10^°(|5m^v|cos2e)^/^ (12) 



Finally, the effective asymmetry has the form: 



a) -T]o , 



(13) 



where r]o = - rj^ - rj^^ - r]^,- 

Note that merely V depends directly on the momentum y. Also note that the 
system of Eqs. (H) is completely equivalent to Eqs. ([|), i.e. we have not used any 
approximations yet. 

3 Evolution of the lepton asymmetry 

We now briefly discuss the different regimes of the evolution of the lepton asymmetry. 
Let us first sketch the reason why the lepton asymmetry can change at all. Due to 
the mixing of with v^, the number densities of the z/g and Vf. can alter. The sterile 
neutrinos that are thus created have no effect on rj. Now if 7^ initially, the 
and Ve number densities will change differently, thus changing rj. Thus the neutrino 
mixing works backreaction on r]. 

In Fig. ID we have plotted an example of the lepton asymmetry evolution as a 
function of the inverse background medium temperature for the case {5m? ,sm29) = 
(— leV^,5 X 10^^) and an initial effective asymmetry 770 = —10^^*^. Figure ^ shows 
the evolution of the resonance momentum ?/res, which is given by 



where C{q) = VQr]{q)q'^^^ / {2Q), and ± is valid for neutrinos and anti-neutrinos, re- 
spectively. For definiteness, we have taken rj to be positive. 

7 




(14) 
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Figure 1: Evolution of the lepton asymmetry. The and T~^^^^ hues were calcu- 
lated numerically with our improved expansion, see section ( ^.2|) , with and without 
the counter-term, respectively. For x > 0.3 we have continued the power-law behavior 
for for T'^ the power law already freezes out at x ~ 0.3. For comparison, we 

have sketched the T^^ power law. 
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Figure 2: Evolution of the resonance momentum. We treat the same cases as in 
Fig. |. 
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/,From Figs. and ^ we can see that there are three different stages of lepton asym- 
metry evolution. During the first stage the asymmetry is so small that the variable 
C{q) can be neglected in Eq. The neutrinos and the anti-neutrinos simultane- 

ously pass through the resonance. Furthermore, any initial asymmetry decreases due 
to the negative back-reaction. 

The second stage is the stage of exponential growth of the lepton asymmetry, 
which occurs when the back-reaction becomes positive. From the point of view of the 
resonance condition, the coefficient C (g) becomes important so that the two resonance 
conditions for u and u are driven in opposite directions (see Fig. One of them 
(in our example u) quickly passes through all the momentum modes to very high 
momenta where the Fermi distribution is negligible, the other one {u in our example) 
decreases to small momenta ?/res < 1- Though the duration of this stage is short in 
time, the asymmetry grows many orders of magnitude and reaches values of order of 
10^^. (For a recent analytical investigation on these two regimes, see |17|.) 

However, we should mention that the behavior of the asymmetry in the exponen- 
tial regime is still unsettled. Several publications |]l8l have claimed that the asymme- 
try starts oscillating in this regime, thereby changing its sign. We believe that this 
effect appears when the momentum distribution is neglected, since then all neutrinos 
pass the resonance simultaneously. Still, there might be a region in the parameter 
space {6m'^,sm26) where these oscillations also occur when the momentum distribu- 
tion is taken into account f^. In the present paper, we will not further investigate 
this question and will assume that no such asymmetry oscillations occur. 

In the third stage the lepton asymmetry grows according to a power law, rj oc T". 
Unfortunately, in the literature there exist different results on the exact power a. The 



two main results are a = —4, see e.g. |]I9|, and a = —1, derived in 0. In the case 
a = —4, the resonance momentum is constant (see Fig. §). Since in this case the at 
the resonance momentum will soon be fully excited, the resonance momentum of the 
u will in reality slowly increase, thus passing through the whole spectrum. Therefore, 
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— 4. In the end, the sector will essentially be excited, so that rj = C(O.l). 
If a = —1, the resonance condition passes much faster through all the momentum 
modes so that the final asymmetry value remains small, r/ ^ 1. In other words, the 
Us are only excited to a fractional amount. 

The rj oc T^^ law arrives due to the counter-term in Eq. (51) of 0. If this counter- 
term is neglected, we find a different power law for the asymmetry, a = —11/3, a 
result which also found when neglecting their counter-term. We also show this 
power law in Fig. |1]. This power law is quite similar to the one with a = —4. 

Evidently the final amplitude of the asymmetry strongly depends on the exact 
power-law dependence. It is therefore very important to find the correct answer. 

4 Discussion of the paper by Dolgov et al. 

We now discuss the treatment of the system of differential equations as done by 
Dolgov et al. in [0. For convenience, we will denote references to equations in their 
paper in the form (D.51). Furthermore, to avoid confusion, we will strictly use their 
notation in this section. 

First we would like to pay attention to some minor misprints, which however do 
not change the results. First, the sign in Eq. (D.20) is wrong. As a consequence, 
the total sign of Eq. (D.51) changes. However, since the initial value of 6o(0) = +1 
instead of —1, the physical content of Eq. (D.51) remains unchanged. 

Next, in Eqs. (D.47) and (D.48), the signs of the first terms, i.e. with the factor 
F, are wrong. 

Also, a correct derivation of Eq. (D.50) yield^ 

where the total sign is the same as in when taking into account the changed sign of 

Z. Furthermore, the second term has changed sign with respect to 0, and there has 

^In terms of our notation, F, 7 and U must be divided by Q. Furthermore, VZ — D corresponds 
to our Vr]/Q. Finally = -f^ + and ^ ^2(^ _ D^a^f + U^{1 + D^/a^f. 
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appeared a new term. Still, the power law holds, since the third term now acts 
as the counter-term. To see this, we take the square brackets in (|l^) at resonance, 
i.e. when = f/^, and use the fact that 7^ <^ D^, which is certainly true after the 
exponential regime. We get [. . .] = [1 + ^(1 — 2)], which is exactly what we get 
from an equivalent treatment of Eq. (D.50). Thus, the counter-term dominates at 
resonance when 7 < \F\/y/8. 

As a consequence of Eq. (20), Eq. (D.51) becomes 



1^ 

Z dq 



-5Bq 



5/3 



2-2 

(T"^cr^ 



1 + 



(72^2 



47' 



(16) 



We also checked Eq. (D.55), and found it to be 



1 - 



A6HI 



;i7) 



i=i,2 - 2CV4g2/3 + (2 

Still, the general physical conclusions do not change, since the missing factor in the 
second term (tl — 1)2 ~ 1 anyway. We will show later on how to derive this type of 
equation, see Eq. (^4|ff). 

20[] . We agree that 



We would also like to comment on the criticism of p| given in 



neglecting the term ^{A — S) in Eq. (D.28) is not allowed when the asymmetry Z is 
essentially damped to zero. However, the effect of this term is merely to slow down 
the decrease of Z. In fact, it is the asymmetry in chemical potential that prevents Z 



from becoming arbitrarily small, see section ^.2| . In any case, during the exponential 
regime, any information on the amplitude of Z before the exponential regime is lost. 
Thus, one can safely neglect the term in the formalism of when one is interested 
in the final value of the asymmetry. 

As a second point, |^D[| criticized that Dolgov et al. neglected the terms H' /Q and 
L' IQ in Eqs. (D.27,D.28). To check the validity of this approximation, we calculated 
these two terms numerically in the formalism used by |Q, and compared the results 
with the other terms in Eqs. (D. 27, D.28). We found that H'/Q and L'/Q were 
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indeed negligible at all times. The differing results in ||20| appear due to the neglect 
of the momentum distribution. We are of the opinion that one cannot neglect the 
momentum distribution, since this would imply that all neutrinos were at resonance 
at the same time, which is simply not the case. 

Third, the author of pO| suggested that the T^^ power law obtained in [^j stemmed 
from neglecting the influence of the lepton asymmetry on the fermion distribution 
function /eq. By expanding /eq in Z ^ ({3)/Tr'^ x lO^^r/j^^, we get 

/eq(Z) ^ /eq(0) I 1 ± " ^2^^^) ^ 



/eqv-y ^^^^"M C(3) X iQio 7r2 l + e-y 

where ± denotes neutrinos and anti-neutrinos, respectively. The correction solely 
leads to an additional term in A' / Q which is given by 

- 2-fZ ^ . (19) 

We compare this with the first term in A'/Q, i.e. —FL, using the leading order in 
Eq. (0): 

lOi^FVt/ , , 

- 20'^ ■ (20) 

One can clearly see that the additional term can safely be neglected near the resonance 
in the power-law regime. 

Finally, we discuss the problems we found in the analysis of . The first occurred 
when the authors diagonalized their set of equations for the symmetric functions, 
Eq. (D.33). When diagonalizing and expanding it in F, they get as the second 
eigenvalue 

"^-"-^■'^ na^-ivzr] • <"> 

where cr^ = 7^ + f/^. Here, the denominator of the second term diverges close to the 
resonance in the case that VZ > 7. We found that this condition is already fulfilled 
before the end of the exponential regime. Furthermore, we found that this divergence 
continued to cause problems in the subsequent analysis. 
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The second problem occurs when F = O {'-/). Then the expansion in F is no longer 
valid. Unfortunately, the counter-term in Eq. (D.50), i.e. our Eq. (|15D , becomes of 
order 1 only when F ^ 7, so that one cannot say whether the counter-term is physical 
or simply an artifact of the expansion. 

Finally, the validity of Eq. (D.49) breaks down when Q'-f <^ 1. Thus, assuming 
that the counter-term is of order 1, (D.49) is invalid for QF <^ 1, which corresponds 
to sin 26 ^J\6m'^\ cos 26 < 1.8 x 10"^ 

5 Improved expansion 

In this section we develop a similar approach to that of using our own notation 
unless quoted explicitly. However, we want to circumvent the divergent behavior 
that occurs in their analysis. To this end, we treat the neutrinos and anti-neutrinos 
separately. 

The system of equations can be rewritten in matrix form as V = A4V: 



( 


\ 




( ° 








^ \ 


('\ 


a' 









-27 





-F 


a 


h' 












-7 


w 


h 


\ I' 


) 




I -F/2 


F/2 


-w 


-7 / 


V I ) 



(22) 



We can apply a non-degenerate transformation 

V = XB . (23) 

For the new vector B we get the equation 

B' = [x-^MX + X-^X') B . (24) 

Now, if A" is a matrix, constructed from the eigenvectors of Ai, then the first term 
on the RHS of Eq. (^) will be diagonal. In the case where is a constant matrix, 
the eigenvectors will also be constants and the second term on the RHS of Eq. ( ^4]) 
will be zero. Then the equations for each component of B will be independent of each 
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other and could easily be solved. Then the general solution of the matrix equation 
will be a linear combination of the eigenvectors of Ai with coefficients B. 

In the general case, when Ai is not constant, there is no straightforward way to 
solve the matrix equation. But in our case in some regimes we can hope that the 
second term on the RHS of Eq. ( p4D will be a small correction in comparison to the 
main term, and then we can solve Eq. perturbatively. 

To begin, we write down the exact eigenvalues of A4: 



-F2 


- 1^2 + ^2 + ^ 


/-4F272 + (F2 + W^ + 72)2 




-W^ + -f^ + ^ 


/-4F272 + (F2 + Vr2+^2)2 




- W^2 + ^2 _ ^ 


/_4F2^2 + (^2 + W^ + 72)2 


-F2 


- Vr2 + ^2 _ 


/-4F272 + (F2 + W^ + 72)2 



= -^-7=2 

/i4 = -7 + ^\/-F2-Vr2 + y _^_4^2y + (^2 + yf/2 + ^2)2 , (25) 

v2 

Note that at resonance W = 0, while out of resonance W ^ -F, 7- We will consider 
three different regimes, F ^ 7, F ~ 7 and F ^ '-f. We will concentrate on the first 
regime and will afterwards stress the problems occurring in the two other regimes. 

5.1 Expansion for <^ 7 

In this case we expand Eqs. (^) in F. We have done so to second order: 

F27 



/ii = -27 + 

fl2 = 



2p2 

F27 



2p2 ' 

/^3,4 = -7T^W^(l + ^) , (26) 

where we have used = 7^ + W"^. 

As a next step, we can construct the eigenvectors, and thus the matrix X, from 
them. Then we re-express the original wave function vector (s, a, h, I) through the 
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new functions bj{T) according to (p3|), again keeping only terms up to order F^: 



2p4 



fei + 62 



4p2 



-— [7 (63 cos — 64 sin Q) — W (63 sin f2 + 64 cos Q)] 



(27) 



h 



F 



— ^ [7 (63 cos — 64 sin f2) + 14^ (63 sin + 64 cos fi)] 
FW 



2p^ 



2p2 
F7 



+ (63 sin + 64 cos fi) 



F7 . , 
61—— - b2—^ + (&3Cosfi - 64 sin n) 



(28) 
(29) 
(30) 



Here, we have at the same time introduced the fast-oscillating variable fl, with 
Q' = W. Also, while the functions bi and 62 are the eigenfunctions of the eigenvalues 
pi and p2, respectively, the functions 63 and 64 are superpositions of the eigenfunc- 
tions to the remaining eigenvalues: the actual eigenfunctions are 62 (cos Q + i sin f2)-|- 
bsl— sin Q + i cosfi) and 62(cos Q — i sin Q) + b^l— sinQ — i cos Q) for p3 and p4, re- 
spectively. We have chosen these definitions to circumvent imaginary values. 



Inserting these definitions into the system of differential equations (^2]) we derive 



b' 



b' 



b' 



2p' 



-bo 



-762 - -bo 



-763 



fFWV fFlV 
1—^1 sinfi + I — 1 cosQ 



f FW 



cos Q 



P 

'F-f' 



sinQ 



p- 



(31) 
(32) 
(33) 



where we have neglected all terms of insignificant order, and 61 can be neglected 
totally. Furthermore, the initial values are given by bo{0) = +1 and 61,2,3(0) = 0. 
Following the corrected Eqs. (D.45-D.51), we get 



/(r) 



F7 



bo 



1 + 



2p2 



1 - 



2f 

p2 



(34) 



Also in this case we find a counter-term which is similar to that derived by 
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To obtain the evolution of the lepton asymmetry, we also have to include the 
anti-neutrinos, i.e. 



2p2 I p2 



Then, by using the differential equation for a and Eq. ([I3|) , we get 

dr] \F\ 



(35) 



Hdy y [Uy,^^)l (f) - /eq(i/, /.)! (f )] . (36) 

•J 



dq 2C(3) 7o 

Note that one power of y cancels when we change from d/dr to d/dq. Furthermore, 
the 7a term in the differential equation for a disappears due to the conservation of 
leptonic charge 

j dyy [feqiy, fJ-ha - /cq(y, /i)7a] = . (37) 

Equation ( p6D is valid in the case F < 7. For -F > 7, the perturbative expansion 
breaks down. We should also check the validity of Eq. (D.49), i.e. the approximation 
of the integral. This approximation is valid when 7 3> 1, which is a weaker condition 



than F < 7 if F > 1, i.e. sin 26'y^|5mgvl cos20 > 1.8 x 10"^. Thus, in the case of 
small mixing angles and small (5m^, the condition that Eq. (|3^) is valid is given by 
7 > 1. From solving Eq. (|36|) numerically, we found that these conditions hold until 
the power-law regime is established for sin 20 < 5 x 10^^ when = — leV^. If we 
neglect the counter-term in Eq. (0), we get the same solution as [|l9l, Eqs. (79-80), 
up to some minor corrections. 



To compare our result Eqs. (pi|-p5|) with 0, we neglect the chemical potential fi, 
use 60 ~ ^0 and expand (/ — I) in Vrj. We get 

. ,38) 

2(t/2-72) F2(5f/4-26f/V + 97^)l\ 



where cj^ = [/^ + 7^. We compare it with the expansion of ([T 

in 2F-yUVZ^ ( F2(f/2-372) , , 

2xlO-^°L ^ (39) 
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This equation is given in terms of the variables as defined in P]. Still we can compare 
the two equations, since the differences to our definitions cancel. We see that, up to 
some minor corrections, the main difference is in a factor of 4 in the F"^ terms. This 
factor has no fundamental consequences. 

5.2 Discussion of <^ 7 

Now we discuss the solutions of Eq. (^) in the region of its validity in more detail. 
First, let us consider the region where the asymmetry decreases (see Fig. |l]). In this 
region, the corrections to the main terms in Eqs. (Q) and (|35D are negligible. We 



expand the integrand of Eq. ( p6| ) in terms of rj and ji/T and get 

feq{y,p)l - feq{y, - p)l ^ fea^iv) ^^^^ ^^^^ " /cq(^) ^^2^~ ( "J + ^ 

P P \P P , 

2 I p2 + p2 



+ Uy)^^^ ( - + - {3[1 - UV)] - 0.2} + 0«, r/.^r^) , (40) 



and b- = bo - bo . (41) 



where we used the notations 

, bo + bo 
^- = ^ 

On the RHS of Eq. (^), the first term corresponds to the main term in Eq. (D.50), 
the second term takes into account the A — S contribution neglected in Eq. (D.50), 
and the last term takes into account the non-zero chemical potential both for the 
7-terms, see Eq. (|), and for /eq, which we rewrite through r]^^ ^ 1.5^. Note that we 
have used fi = — /i, which is valid for large temperatures. Before substituting Eq. ( ^OD 
into Eq. (^), let us write the differential equations for b±, in which we can neglect 
the small difference between 7 and 7: 

db+ _ _F^ ( -r^ + + {Vr]f UVt] \ 

^ . 2!±4±(M,_ + 1^,,) . ,42) 

dr 2 y p^p^ p^p^ J 

In the case of small asymmetry rj we approximately get from these equations 1 
and 6_ oc sin^ 291]. Using this fact, we can substitute Eq. (^) into Eq. (|36D and get 
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the following equation: 

^ = -Ct^ [h{q) - sin' 26 h{q)] + r^uJM , (43) 

where the functions Ii{q) are integrals over the momentum which do not depend on 
the asymmetries. The term Ii{q) corresponds to the main term in (D.51), and without 
I2 and I3, the asymmetry decreases to infinitely small values. 

The second term has opposite sign and is proportional to an additional power of 
sin^ 26. For sin^ 26 ~ 10~^, this term becomes of the order of the first one, but it can 
not cancel it completely. We would like to mention that, although this term seems 
to be related to the expansion in F, it is actually of a different nature, and should 
therefore not be neglected here. This term corresponds to the A — S contribution in 
terms of and it is proportional to h_ in our notations. 

But the actual reason that the lepton asymmetry does not decrease to very small 
values is the non-zero chemical potential, which corresponds to I^. When the asym- 
metry, which evolves according to Eq. (|^), becomes small enough, 77 = 2riij^ — rj^ <^ 
f]uc ~ ~''7o/2, the last term in Eq. (|^) starts to compensate the first two terms, and 
prevents the asymmetry to get tiny values. 

We will not consider the regime of exponential growth of the asymmetry (see 
Fig. |l]) separately, let us just note again that during this stage the resonance condition 
for anti-particles (particles) go to very large momenta, while the resonance condition 
for particles (anti-particles) go to small momenta, ?/ < 1, for 77 > (r^ < 0), see 
Fig. ^ Therefore, in the power-law regime, one may neglect either neutrinos or 
anti- neutrinos in (|36D, since the contribution of one of them is negligible due to the 
exponential damping of the Fermi distribution function. 

Finally, we discuss the power-law regime. In this regime, we can neglect either 
neutrinos or anti-neutrinos. For definiteness, we consider the case 77 < 0, where the 
resonance momentum of the anti-neutrinos is very large. Thus we can neglect / in 
Eq. (^). Also, we will assume that the asymmetry is still 77 -C 1, so that we can also 
neglect the chemical potential. This is in contrast to the case described by Eq. 



where the main terms of the v and v almost canceled each other, so that the even 
the small chemical potential would become important. Here, however, the v have no 
influence at all, so that no cancellation of the main terms appears. Furthermore, we 
assume 60 — I5 which is valid for small -F/7. To be able to integrate analytically, we 
expand W around j/res, given in Eq. (0): 

W^(l/)~CH/(l/-2/res) + ..., With = + V^r/ . (44) 



Then the leading order in (^6]) gives 



dr] f J r / \ 



dyfeqiy)- 



dq 4C(3)g2y + 

Jeqiyres)-^ / dy- 



4C(3)g2 ny,,,)_l 1 + _ ^^^^2 



^rcs 

00 

res I 



"^'^/eq(2/res)^#^^ / dz-^^ , (45) 



4C(3)g2-'^-'^-y(2;,,3) cw J l + z 

— oo 

where we have substituted z = cwiy — |/res)/7res- The final integral gives vr. Analo- 
gously, we find the next order corrections, so that we get 

dq ACiSr^'^y'''^ V Snyrcs)) ' ^ ^ 

which is our differential equation for the asymmetry rj in the power-law regime. We 



can compare this to Eq. (|T^, which is the corresponding equation for the approach 
by Dolgov et ai. The interesting point is that the counter-term, compared with the 
leading term, is the same in both approaches up to a factor of 2 (the counter-term in 
Eq. ( p!?!) corresponds to F^/(167^(?/res)) in our notation). 

/^From Eq. ( ^61 ) it is easy to derive the power law for rj. Let us first neglect the 
counter-term. If we further assume that C{q) ^ q in (p!4D, we get y^es ~ <5/ iVoVQ~'^^^) 
and can neglect the first term in Cw Then 



dr] -F'tt ^ ^ O 



Uyr.s)^,r\"'' ^ C,r\"' ■ (47) 



dq 4C(3)''^'^'^'=''V^o' 
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If we multiply by dq and integrate, we get 

V=[yiC,) g"/^ (48) 

Here, we have neglected the boundary conditions, which become negligible for large 
q in any case. We see that this corresponds to the 7] oc T^^^^^ power law found in [Q 
when neglecting the counter-term. We should however mention that even for all the 
simplifications made above, this power law is only approximate, since contains the 
Fermi distribution feqiVres)- Since y^es oc g^/^^ for the derived power law, /eq(?/rcs) will 
slowly decrease, thus changing Eq. (^). 

If we take the counter-term into account, the condition for the power law is that 
the brackets in Eq. ( ^61) equal 0. Using again the approximation y^cs ~ <5/(Vb'7Q'~^''^), 
we can derive 

Q fss^V" ,,4 

which is exactly the rj oc T^^ power law found in 

Let us still have a look on the evolution of / for the 1] oc T^^ power law. In Fig. ^ 
we have plotted I as given in Eq. (0), once with and once without the counter-term, 
for the case where F = a/27 ^it resonance. Then the counter-term merely compensates 
the leading order term near the resonance, but the integration over the momentum 
still yields a similarly large value when taking the counter-term into account. Thus, 
it is safe to say that the counter-term is of minor importance for F < v^7- This 
corresponds to the result in Eq. (^6|). 

To summarize, this ansatz gives similar results as the approach of 0, so that 
the divergence in their approach seems not to be critical. However, our approach 
has shown two further advantages: First, we need not make the adiabatic approach 
which they used for L and H. Second, our result yields the correct behavior for very 
small rj before the exponential regime, since we included the term {F/2){a — s) and 
the chemical potential. 
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Figure 3: l{y) for the case F = V27 at resonance and rj oc T ^. The abscissa 
shows the momentum relative to the resonance momentum, i.e. y/Ures- The dashed 
hne shows the first term in Eq. (0), while the solid line shows the full equation. 

5.3 F>7 

This region can occur in two cases: either when F, and thus the mixing angle sin 26', 
is large, or when 7 is very small, which is the case for r ^ 1. Since we only discuss 
small mixing angles in our analysis, the first case is not important. However, the 
second case appears, since for large asymmetries, the (anti-)neutrinos at resonance 
have r = g/y, where g > 1 is large and ?/ ~ 0.1 is small. 

In this limit the eigenvalues Eqs. (|25|) should be expanded in 7. However, in this 
regime, the approximation in Eq. (D.49) does not hold, since 7 is no longer very large. 
Therefore, we did not investigate further in this case. 



5.4 



7 



In this case let us rewrite 7 = F{l + x) and suppose that x ^ 1. Then we can expand 
the eigenvalues (|25|) in x- To the main order, which does not contain x, we get: 



-F ± -^^-W^ ± WVW^+AF^ + 0{x) ■ 
v2 



(50) 
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The important feature of Eq. (|50D is the fact that at resonance, W = 0, all 4 eigen- 
values coincide. This means that the resonance matrix Ai changes its rank, and thus 
diagonalization becomes impossible. 

We conclude that the approach of Dolgov et al. is useful and correct as long as 
F < 7, but that it seems very difficult to develop this kind of approach for F ^ 7. 

6 Solution for |a — s — 1| ^ 1 
6.1 Solution for one momentum mode 

In this section, we present a different analytic solution of the system of differential 
equations (|) in the power-law regime for a given momentum y. We introduce two 
simplifications. First, to circumvent the non-linearity, we give the lepton asymmetry 
an explicit form, 

r^{T) = r^,T^ , (51) 

where a represents the power-law behavior. Second, we treat the change of the 
variables a and s, expressed by the function /(a_s) = (a — s — 1), as a perturbation in 
the differential equation of /. In terms of the physical content, f[a-s) ^ 1 means that 
there is only very little conversion from active to sterile neutrinos. The advantage of 
this approximation is simple: provided that the power law holds, the condition 
f{a-s) ^ 1 holds for all timesQ. 

To start with, we write down the differential equations for h and /: 

/' = ^{l + ha-s))-Wh-^l , (52) 
h' = Wl-^h . (53) 

Subsequently, we will assume that f[a-s) is an explicit function. Then the above 
equations form a complete set of differential equations for h and /. 

^This condition corresponds to — 1 <C 1 in the formalism used by and others. 
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By solving (|52D for /i(r), taking its derivative h'{T), and inserting both into (|53|) , 
we find the differential equation of second order for /(r): 

Wl" + {2W-f - W) I' + (W^ + -f^W + W-i' - 7iy') / 

^ '(1 + (l^'-7W^)-H^/(a» J = 0. (54) 



+ 2 



Analogously, we find 



Wh!' + (2Vr7 - W) h' + {W^^ + 7^1^ + W^y - -fW) h 

-l{l + f^,_,^)W' = 0, (55) 

where the only difference to ( ^If ) is given in the inhomogeneous F-term. 

We can now solve the equation for / analytically. To this end, we first write down 
the solution for the homogeneous differential equation, i.e. F = 0: 

T T 

IhoJr) = I ki cos( / W2dT2) + k[ sin( / W^s^ts) 



X exp y- j 72(ir2 j . (56) 

Here, the coefficients k\ and A;^ are given by the boundary conditions at time r^, and 
the subscript 2 means that the parameters are evaluated at T2. The solution for the 
inhomogeneous differential equation can then easily be found, e.g. by using [PT|] , 

F } (1 + /(a^.))i jWj - 7iiyi) - W^i(/(;_.))i 



Anhom(^) 



X sin I - y W2dT2 I exp I - y 72^X2 j dn . (57) 



The complete solution is then /analytic = Ihom + Anhom- 

Analogously, we get h = /ihom + ^iinhom, where /ihom is identical to the homogeneous 
solution for /, (^), now with coefficients k^ and k^. The inhomogeneous solution is 



h -I 

"■inhom 



+ f(a-s))i sin - y W2dT2 exp - y 72^X2 dn . (58) 
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Now we turn to the variables a and s. It is straightforward to solve the differential 
equations describing them, see (||), provided that / is an explicit function: 



a r 



; exp I - y 72(ir2 \ - F j ^expi(Ti) exp I - y 72(^X2 j dri (59) 



s{t) = Si + Fj /expi(ri)c/ri . (60) 

We still need to settle the boundary conditions: From |0 we know that for Tj = 0, 
h = I = a = and s = — 1. It follows that k\ = k\ = = k!l = ai = Q and Sj = — 1 
as initial conditions. 

The next step is to connect the two systems of differential equations, (/i, /) and 
(a,s), using perturbation theory. The leading order, labeled with the subscript 0, is 
defined by /(a-s) = 4xpi = 0. We find a^ij) = so(t) = 0, while /q and ho are given by 
(13) and (H), respectively, using f(^a~s) = f[a-s) = 0- 

The first-order corrections, Oi and si, are found by inserting /gxpi = lo into 
the analytic solutions, ( pOj ) and (|60D, respectively, and for h and / by inserting 
/(a-s) = cii — Si — 1 into ( ^8|) and (|57|), respectively. For higher order corrections, this 
procedure is repeated recursively. It is clear that this expansion holds only provided 
that /(a-s) ^ 1 for all earlier times. 

6.2 How to check the consistency of a power law 

In the previous subsection, we derived an analytic solution for the density matrix of 
a single momentum y. The analytical equation for rj in the power-law region is found 
by using the differential equation for a and Eq. (p!3D, 



where we have neglected the anti-neutrinos and have used the conservation of leptonic 



charge, Eq. (|37|) . We can insert the solution of / as given in the previous subsection 
by using the relation t = q/y. However, note that the resulting equation for rj is valid 
only as long as /(a-s) <^ 1- 
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Furthermore, we need to define W, which is the only parameter in the analytic 
solution of / that is dependent on r] and y. If we define 77(g) = rjoq'""^^, where a is 
the temperature behavior of rj, rj cc T", we can write 

W{t, y) = U{t) - Vor7or-(^+")/=^y-(^+")/=^ • (62) 

Having given W in this form, we now have two unknown parameters: the power-law 
behavior a and the asymmetry factor r^o- 

Checking the consistency of a power law is now simple in principle: given the 
correct values for a and r/o; the integration of the rhs of ( pT| ) for any q in the power- 
law regime should yield 



-? m q''-'^" . (63) 



, ^0. ) analytic ^ 

6.3 Results 

We now consider the case derived by i.e. a = —1. In this case, ( p^ ) does not 
depend on y directly, while the y-dependence persists indirectly through r = q/y. 
The remarkable consequence is that Tj-^s is independent on y, and thus is the ratio of 
F to 7 at resonance, R = -F/7res- 

We found it convenient to substitute the integration over y in Eq. (0) with an 
integration over r: 

dr] Fq^ f dr 



M 2C(3) y r3-^-'«/^"M ■ '«^) 

Also, we use the parameter R instead of 779. The advantage of this is obvious: the 
power law in [0] originates in the counter-term compensating the leading order term. 



which occurred when R = const = 0{1). This yields Tres = yF/{R6). Furthermore, 
by solving W (Tres) = 0, we get 

\r]o{R)\ = QV,-\UR) [r'HR) - 1] • (65) 

The following calculations were done using the fits presented in Appendix. To 
begin, we have calculated the evolution of /analytic using f(^a-s) = for 0.1 < -R < 1000. 
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We have given an example for R = 100 in Fig. Within this range, / showed a similar 
behavior, the only differences occurring in the scaling of r and the amplitude of /. 
Then we calculated the first-order contribution from f{a-s)- In terms of the formalism 
used in 0, this corresponds to including contributions from 6q, i.e. the origin of the 
counter-term they found. We found that the correction was negligible for R < 10, 
and that it became of order of the main term only when R ^ 1000, and also then 
stayed smaller than the LO term. 

Thus, we conclude that the counter-term found in |^ must be an artifact of the 
expansion in F, which is only valid for R <^1. 

Looking at Fig. ^ one might argue that even without the NLO term, the oscilla- 
tions after the resonance might be the origin of a large counter-term. However, we 
found this not to be the case. We calculated Eq. (|64D for 0.2 < i? < 10. According 
to (|63D, the result should be equal to r]oq~'^^^/3 for all q, with t]o as given in (p5|). 
We found that this was not the case, even when including the first-order contribution 
from /(a-s)'- as an example, we show r]oq^'^/^ /3 and ( |64D in Fig. |, with R = 1. 

By the way, we found that ( ^ ) could be fitted very well by the simple 

V "<? / analytic 

function 



. , , 5.4 X 10-^ ye J ^ , (66) 

V^^/ analytic ^ ^^-s/ 

in the region 0.2 < i? < 10, 0.1 < q < 60. This corresponds to our understanding 
that only the momenta close to resonance are significant. From this fit, we see that 
( ^) has a totally different q behavior than the expected 

X'^Q / analytic 

The reason that the oscillations do not bring a counter-term can be seen in the 
evolution of s, shown in Fig. |^. The better part of the asymmetry is created in the 
region around the resonance, where / varies slowly. The oscillations after this region 
very soon become fast, and each half period approximately compensates the previous 
one. 

For the case a ^ —1, the situation is more complicated, since then the evolution 
of / depends on y. We intend to eventually present the results for this case in a later 
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Figure 4: /(r) for the case R = 100, a = —1. We have plotted the numerical solution 
(solid line), the analytic solution with f(^a-s) = 1 (dashed), and the analytic solution 
with first-order correction (dashed-dotted). 
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Figure 5: s(r) for the case R = 100, a = —1. We have plotted the numerical solution 
(solid line) and the analytic solution with f^a-s) = 1 (dashed). 
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Figure 6: (^^) ^ . (solid) and the expected value rjoq ^/^/3 (dashed) for the case 



^ d'Q I analytic 

= 1, a = -1. 



publication. 

However, we decided to consider here the special case a = —11/3, i.e. the power 
law derived in Eq. (|48|) , numerically. To this end, we assumed this power law, with 
Vo = (9/llC^)^^'^, where was defined in Eq. (|^). Then we calculated the evolution 
of / for a grid of momenta y numerically, using the exact differential equations (P), 
and integrated the results for several q according to Eq. (|6TD . Then we compared the 
result with the expected value, in this case given by 

fdr]\ 11 



dq 



9 



■Voq 



2/9 



(67) 



expected 

We have plotted the ratio in Fig. ^ If our result was perfect, the ratio would have 
been 1 for all q. The discrepancy has several origins: the result depends on a and 
rjQ. As an example, choosing a different rjQ will shift the curve upward or downward. 
Hereby, already a change of a few percent would suffice to make the curve in Fig. |^ 
cross unity. Furthermore, a change in a alters the curvature of the plot. In particular, 
the downward slope at small q is due to the fact that we have used the approximation 
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Figure 7: (^) /(t^) as a function of q for a = —11/3. For simplicity, we 

V "9 / num V "9 / expected 

have set feqiVres) = 1/2) which is a good approximation for not too large q. 

C{q) ^ q when deriving the power law, see Eqs. (^6|-^8D. In this context, the result 
presented in Fig. ^ is fairly good. We conclude that it is safe to say that the real 
power law is around a ~ —4. 

7 Conclusions 

In this work we have analytically investigated the system of kinetic equations which 
governs the creation of lepton asymmetry in the Early Universe due to active-sterile 
neutrino oscillations with small vacuum mixing angle sin2^^ < 10~^ and negative 6m^. 
We have improved the analytical approach of 0, taking into account the non-zero 
chemical potential of neutrinos and relaxing the adiabatic approximations made in 
i- 

In analogy to [Q, we have derived an equation for the asymmetry evolution 
[Eq. (|36D], which we found to be valid only if the collision terms dominate the ki- 
netic equations (H), or in our notation if F < 7. Except for very small values of the 
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asymmetry, our approach gave the same results as found in Also, the main term 



in Eq. ( p6|) coincides with the so-called static limit derived in [jT9| . However, Eq. ( pq) 
contains an additional term which compensates the main term when F ~ 7. 

It was found by the authors of that this additional term resulted in a growth of 
the asymmetry according to a power law. However, we found that the additional 
term in Eq. ( ^6|) becomes important only when the expansion breaks down, or F/j = 
(9(1). But strictly speaking, this fact does not prove that rj oc is wrong, because 
the additional term already exists when the perturbative expansion is valid. 

In order to check the validity of the T^^ power law, we developed a new analytical 
approach which is independent from perturbations in -F/7. We expand the kinetic 
equations not in the coefficients (as done in but by expanding in the combination 
of functions (a — s — 1), which remains small as long as the number density of the 
sterile neutrinos remains small for all momenta. 

This condition strongly depends on the power-law index. For the T^^ power law, 
the number density of the sterile neutrinos remains small for all momenta. This fact 
allowed us to check the consistency of r/ oc T^^, and we found that this power law is 
not a solution of the kinetic equations. Thus, we proved analytically that rj cx is 
erroneous. 

For the power law, which corresponds to a full transition of active to sterile 
neutrinos, the condition (a — s) ~ 1 breaks down already in the region F/7 < 1. For 
the T~^^^^ power law, this condition holds until F/7 ~ 1. Thus, our new method is 
not useful beyond this point. 

We conclude that our Eq. ( p6D provides a correct analytical solution, which is in 
qualitative agreement with the numerical results of 0, ^ |^ in the region prior to the 
formation of a power law. For sin 20 < 5 x 10~^, our solution shows the formation 
of the power law T^^^^^. On the other hand, we can not apply our Eq. ( ^6|) in the 
region where the MSW effect dominates and can not decide if this power law changes 
to or not. Either way, the final value of asymmetry should be large. 
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When this work was nearly complete we received a manucript from P.Di Bari, 
R. Foot, R.R. Volkas and Y.Y.Y. Wong who also conclude that the power law 



1] (X T ^ is erroneous by showing that the approach of does not work in the region 
where the MSW effect dominates over the second-order effects. 
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Appendix: Fits for / 

Since the exact analytic function for / is an integral, and we are forced to integrate 
over /, we have fitted / to speed up the numerics. Here we present the fits used in our 
calculations in subsection |6.3| . 

At very small r, the adiabatic solution, given by h' = V = 0, is a good approx- 
imation. However, for small 7, the condition /I'adiab W^^adiab, 7^adiab breaks down 
long before the resonance. We therefore relax the condition h' = 0, while we maintain 
/' = (which is valid much longer due to the term F/2). Thus, by setting /" = /' = 
in (|5^) , we find 

_ F ^W-W 

The error of this solution is of order /semiad/(-^/2)- 

For the region before and during the resonance, we decided to calculate the ana- 
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lytic solution of / for a sample of points and fit a function to it. We found that before 
the resonance, / showed an exponential-like behavior, which is why in this region we 
used a fit of the form 

^expfit = e°'^'^\ where Q;(r) = + aix + . . . (69) 

Using n = 30 points and a polynomial of order (n — 1), we found it necessary to do 
two such fits to achieve an accuracy of better than 1%. 

Around the resonance and shortly after, we applied two linear fits, also with 
n = 30. 

Finally, we fitted the oscillating behavior after the resonance. We found that it 
is a good approximation to add the homogeneous solution (|56| ) to the semi-adiabatic 

^oscfit(T) = ^hom(T) + 4cmiad(T) • (70) 

Then the semi-adiabatic fit accounts for the inhomogeneous term 

When using the analytic solution for r > tq, we found that the integrand diverges 
at the resonance, ti = tq. However, this problem can easily be circumvented: between 
Tq and some point ^ tq we used a quadratic fit for / with boundary conditions given 
by the analytic solution at Tj = tq, while for r > we used the analytic solution with 
boundary conditions given by the quadratic fit at Tj = r^. We took = r[W^(r) = F]. 

We should also mention that for r ^ tq the integrand of the analytic solution is 
oscillating fast, so that it is difficult to calculate the integral starting from r = 0. 
Therefore, we have used the semi-adiabatic solution up to some Tj <^ tq, where the 
semi-adiabatic solution still is a very good fit, e.g. has an error of less than 10~'^. 
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